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We report a study of the complex AC impedance of CVD grown graphene. We mea- 
sure the explicit frequency dependence of the complex impedance and conductance 
over the microwave and terahertz range of frequencies using our recently developed 
broadband microwave Corbino and time domain terahertz spectrometers (TDTS). 
We demonstrate how one may resolve a number of technical difficulties in measur- 
ing the intrinsic impedance of the graphene layer that this frequency range presents, 
such as distinguishing contributions to the impedance from the substrate. From our 
microwave measurements, the AC impedance has little dependance on temperature 
and frequency down to liquid helium temperatures. The small contribution to the 
imaginary impedance comes from either a remaining residual contribution from the 
substrate or a small deviation of the conductance from the Drude form. 

PACS numbers: 78.67. Wj, 78.70. Gq, 07.57.Pt 



I. INTRODUCTION 

Graphene is a material consisting of a sin- 
gle atomic carbon layer arranged in a 2D 
honeycomb lattice, discovered in 196£ffland 
studied extensively since therPby the surface 
science community. It has attracted wide- 
spread interest for both its novel electronic 
properties and Dirac band dispersion as well 
as its broad application potentiaP^. Due to 
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its high mobility, it also has been proposed to 
show great promise for high speed switching 
in microwave and terahertz devices^^l anc i 
terahertz plasmon amplification^. 

Recently it has been demonstrated that 
large-area monolayer graphene films can be 
grown by chemical vapor deposition (CVD) 
on copper foils^, following the precipitation- 
based growth of somewhat non-uniform few- 



layer graphene films on Ni foil d 15 * 16 *. This 
method^^ allows the growth of large scale 
graphene films that can be transferred to var- 
ious substrates, which will be essential for 
any practical device applications. The avail- 
ability of large area uniform graphene also al- 
lows access to these materials from a greater 
variety of experimental techniques including 
studies of their long wavelength electromag- 
netic response. Their complex microwave 
and terahertz frequency dependent response 
are of particular interest and their under- 
standing is crucial in order to use graphene 
for fast electronic devices. However, it has 
traditionally been difficult to get significant 
broadband spectral information in these fre- 
quency ranges, particularly in the microwave 
regime. Microwave experimental techniques 
are typically very narrow band and may at 
best allow the characterization of materials 
at only a few discrete frequencies. 

In this study, we make use of our newly- 
developed microwave "Corbino" spectrome- 
ter to measure the broadband microwave re- 
sponse in the frequency range from 100 MHz 
to 16 GHz of CVD grown graphene at tem- 
peratures down to 330 mK. This technique 
allows one to gain broadband spectral in- 
formation in the microwave regime. Mi- 
crowave techniques are typically very narrow 
band. We present data for both the sheet 
impedance and complex conductance. The 



measurement of the intrinsic impedance of 
a single atomic layer film on an insulating 
substrate presents a number of experimen- 
tal challenges in the microwave regime. As 
a non-resonant technique, the Corbino spec- 
trometer requires an intricate calibration pro- 
cedure. More importantly, any attempt to 
measure the intrinsic impedance of graphene 
will be affected by capacitive coupling to its 
dielectric substrate. At microwave frequen- 
cies, the impedance associated with the sub- 
strate can be a substantial fraction of the 
graphene impedance. We detail the manner 
in which substrate effects may be calibrated 
for, as well number of other difficulties 
peculiar to this frequency range that must be 
overcome. We also compare our data to that 
we measure at higher frequencies using time- 
domain terahertz spectroscopy. Both mea- 
surements techniques are capable of measur- 
ing the complex optical response functions 
SIS db function of temperature and frequency, 
without resorting to Kramers- Kronig trans- 
forms. 

CVD grown Graphene was prepared by 
methane and hydrogen at pressures of 1.5 
mbar over a 25 /im thick Cu foil. The 
graphene films are coated with PMMA and 
then the Cu foils are dissolved in an aque- 
ous solution of FeCl3. The graphene is rinsed 
several times with de-ionized water and can 
then be scooped out of solution onto a 380 
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band peaks in their 532 nm Raman spectra 2 ' 
(Fig- 0(a)) 



2000 2500 

Raman shift (cm" 1 ) 



2200 

CO 

E 

.c 

O 2000 
□ 

or 



— R n before annealing 

— R n after annealing for 24 hours 
+ R n after annealing for 24 hours at 260MHz 




100 150 200 250 

Temperature (K) 



FIG. 1. (Color Online) (a) Raman spectra 
of graphene on Si. The Raman signal from a 
bare Si substrate has been subtracted, (b) Re- 
sistance per square Rn for one of the CVD grown 
graphene samples as a function of temperature. 
The black crosses are the microwave data at 260 
MHz. After the same sample was annealed at 
about 60 C in vacuum for about 24 hours, the 
temperature dependence of its resistance shifted 
as shown by the red curve. 



/iin thick clean high resistivity Si substrate. 
High purity Si was used as it has a purely di- 
electric contribution to the impedance of the 
graphene-silicon multilayer. The sample is 
allowed to dry and adhere to the Si and then 
the PMMA is removed by acetone. The re- 
sulting films are verified to be of high-quality, 
predominantly single layer graphene from the 
intensities and positions of the G- and 2D- 



II. MICROWAVE 
MEASUREMENTS 

The microwave 'Corbino' technique has 
been developed to provide the frequency 
dependent response in this spectral range, 
where broadband measurements have been 
traditionally been very challenging. It has 
been used in recent years to give important 
information about high temperature super- 
conductors, heavy-fermion, electron glasses, 
and thin film superconductors^" In this 
technique one measures the complex reflec- 
tion coefficients S™ from a thin film sample 
that terminates a coaxial transmission line. 
Typically, electrical contact is made by press- 
ing a modified microwave adapter against the 
sample with donut shaped gold pads, which 
matches the radial dimensions of the coaxial 
line. The technique necessitates the use of 
three calibration samples with known reflec- 
tion coefficients (20nm NiCr on Si, a blank 
Si substrate and a bulk copper sample) to 
remove the effects of extraneous reflections, 
phase shifts and losses in the coaxial cables^. 
Using the telegraphers' equation, the actual 
reflection coefficient at the sample surface S^i 
is expressed as: 



ous independent study of thin superconduct- 
ing films on identical Si substrates 7 . In that 
study, we used amorphous metal films with 



S?i = 11 . r. (I) ing films on identical Si substrates 7 . In that 
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Here, the complex error coefficients Ed, E$ 

, ,-, 1-1 j j. .. .. scattering rates (~ 100 THz) so high that 

and Er represent the so-called directivity, 

, , , a ,. , . the intrinsic AC impedance of the film itself 

source match and reflection tracking errors 

, rp, , , , £ in the normal state was purely real and could 

respectively, lhey are temperature and fre- 

, , , rr. . . m be deduced from the DC resistance exactly, 

quency dependent complex coefficients, lo 

ii^- i ryeff Thus, Zo ub can be calculated by comparing 
extract out the sample sheet impedance Z s , 

. , , , ,. the measured AC impedance of the metal film 

m principle the standard equation 
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with its known value. Knowing Z§ ub one can 



rreff f + S'ii fn\ isolate the intrinsic broadband impedance of 

Z s = 9 1 _ S a Z ° ( 2 ) 

11 the graphene layer. 

may be used. Here Zq =50 ohms is the 

In addition to measurements at microwave 

characteristic cable impedance and g = 

frequencies, the Corbino spectrometer system 

2tt/ ln(r2/ri) is a geometric factor where r<i 

can measure the two contact DC resistance 

and T\ are the outer and inner radii of the 

simultaneously using a lock-in amplifier and 

donut shaped sample. However additional 

a bias tee. Multiplying the measured resis- 

complications are presented in the case of 

tance between the inner and outer conductors 

materials like graphene where the impedance 

of the coaxial cable by the geometric factor g, 

of the film is comparable to that of the sub- 

we obtain resistance per square. As shown in 

strate, as the substrate contribution must be 

Fig. [11(b), the resistance per square as a func- 

taken into account explicitly. In the thin 

tion of temperature is approximately tem- 

film approximation and under the assump- 

perature independent (~ 3% over the range) 

tion that only TEM waves propagate in the 

with only an upturn below 30K as the princi- 

transmission lines, the effective impedance 

pal distinguishing feature. The sample prop- 

for a thin film of impedance Zs backed by a 

erties are changed only slightly by anneal- 

substrate with characteristic impedance Zg 

_ ing in vacuum for about 24 hours at 60 C 

can be shown to be^i 

degree. As shown in Fig. [ijb) the resis- 

„ ff Zq tance per square increased by 5 % with al- 

^ =7vV- ( 3 ) 

-L + most the same temperature dependence. We 

To isolate the impedance of the graphene ascribe the change in the overall scale due 

layer, we use Eq. K^with Zg ub from our previ- to a change in carrier density by driving off 
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FIG. 2. (Color Online) (a) Calibrated 
impedance with (black) or without (red) cor- 
rection for the substrate contribution at room 
temperature. Corrected (b) impedance and (c) 
conductance (normalized by Gq = ne 2 /2h) as a 
function of frequency in the range of 100 MHz 
to 16 GHz at different temperatures. Different 
colors in both panels indicate different temper- 
atures. In all the panels, solid lines are the real 
parts while dashed lines are for the imaginary 
parts. 

absorbed gases. 

In Fig. [2j we present the results of our 
broadband microwave measurements on one 
particular graphene sample from 100 MHz to 
16 GHz at temperatures down to 330 mK. 
The small oscillations are the residual effects 
of standing wave resonances in the transmis- 
sion line that have been imperfectly removed 



by the calibration procedure. In Fig. [2^a), 
we compare the effective impedance at the 
sample surface calculated from Eq. [2] with 
the impedance of the sample after the sub- 
strate correction described above. One can 
see that the correction becomes significant at 
higher frequency where the effective capaci- 
tance of the dielectric substrate plays a larger 
role. After correction the impedance becomes 
primarily real and frequency independent as 
expected for a conductor with a scattering 
rate larger than the measurement range. One 
can see that it is essential to perform such 
a correction to quantify the impedance cor- 
rectly. 

In Fig. [5] (b), we plot both real [Z\) and 
imaginary (Z 2 ) sample impedance corrected 
for the substrate contribution SIS db function 
of frequency at 5 different temperatures. We 
can see that the frequency dependance of 
impedance at the base temperature of 330 
mK and room temperature are almost the 
same. The vertical difference in those tem- 
peratures matches with the difference of mea- 
sured DC resistivity. As clearly seen from 
Fig. [2] (b), the real and imaginary parts 
of impedance have little dependence on fre- 
quency down to low temperatures. This can 
also be seen in Fig. [I] (b), where the resis- 
tance at the low frequency of 260 MHz at 4 
different temperatures is also plotted. These 
data follow the DC values closely indicating 



a consistency of AC and DC measurements. 
Also, the real and imaginary parts of con- 
ductance have little dependance on tempera- 
ture which means that the Drude response of 
electrons in graphene does not change a lot 
over this wide range of temperatures (from 
room temperature to 330 mK). It also indi- 
cates that the scattering rate r - the average 
time between scattering events - bears little 
dependance in temperature. 
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FIG. 3. (Color Online) Ratio (red) of real and 
imaginary parts of conductance as a function of 
frequency in the range from 700 MHz to 9 GHz 
at 15 K before annealing. The black curve is a 
linear fit with zero y axis intercept. 

In Fig. [2] (c), we plot the complex con- 
ductance obtained by using the inverse of the 
data Fig. [2] (b). Since in the thin film limit, 
the complex conductance is the reciprocal of 
complex impedance, it also has almost no de- 
pendance on frequency. Here we have ra- 
tioed this data to the quantum of conduc- 
tance G Q = ixe 2 /2h = 1 / 16433 Ohms" 1 ex- 



pected for a graphene sample with its chem- 
ical potential tuned to the Dirac point. Its 
large dimensionless scale shows that the car- 
rier density for this sample is high with the 
chemical potential far from the Dirac point. 

Although the imaginary part of the con- 
ductivity is very small in Fig. 2, it is not zero, 
which gives a measurable r from the data. 
Within the Drude model, the complex con- 
ductivity from a charge responding to a time 
varying external oscillating electromagnetic 
field with frequency u is a(u) = t — 
where n is the electron density. Inspection 
of this equation shows that the ratio of the 
imaginary to real parts of 02/01 gives ur as a 
function of u such that r can be determined 
from its slope. In Fig. [3j we plot this ratio 
for the conductances (G = ad) as a function 
of frequency at 15 K. GijG\ at other temper- 
atures give similar results as both G2 and G\ 
have little dependence on temperature. The 
slope of G 2 /Gi gives us an estimate of the 
relaxation time, which is about 25.7 ps. A 
rough estimation of scattering rate is then 
T = 1/r = 38.9 GHz. This is quite close to 
the value for the Drude scattering rate of 36.4 
GHz obtained independently through fitting 
our data using a Drude-Lorentz modeP^. Us- 
ing this extracted scattering rate, we can es- 
timate the mean free path I = Vp * r = 28.3 
//m with Fermi velocity vf = 1.1 * 10 6 m/s ? . 
However, we should note that this small value 



of the scattering rate is interesting as it is at 
odds with that inferred from previous stud- 
ies using higher frequency time-domain tera- 
hertz spectroscopy®' or far-infrared reflectiv- 
ity measurements^. 
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FIG. 4. (Color Online) Real (solid) and imag- 
inary (dashed) parts of normalized conductance 
as a function of frequency in the range from 150 
GHz to 1.0 THz at different temperatures for a 
similar graphene sample. 



III. TERAHERTZ 
MEASUREMENTS 

The small value of the scattering rate is 
also different from that inferred from our 
own time-domain terahertz spectropscopy 
(TDTS) measurements on another similarly 
prepared sample. In TDTS an ultrafast 
laser pulse excites a semiconductor switch, 
which generates an almost single cycle pulse 
with frequencies in the terahertz range. The 
transmitted terahertz pulse's electric field is 
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mapped out function of time. The ra- 
tio of the Fourier transform of the transmis- 
sion through the sample to that of a reference 
(usually the substrate on which the sample 
is deposited) gives the complex transmission 
function of the film under study. This can be 
inverted using standard formulas in the 'thin 
film approximation' f(u) = 1+n +^ {u)d e^ B 
to get the complex conductivity. Here $ s is 
the phase accumulated from the small differ- 
ence in thickness between the sample and ref- 
erence substrates and n is the substrate in- 
dex of refraction. We have measured in ter- 
ahertz frequency range from 150 GHz to 1.0 
THz at temperatures down to 6K with flow- 
ing He4 gas. It is possible that the He4 gas 
environment for graphene sample in TDTS 
might make a slight difference for its scatter- 
ing rate compared with the Corbino system 
where the sample is sealed in high vacuum. 

In Fig. [4j we show terahertz complex 
conductance data for a sample from a dif- 
ferent batch. This sample has a higher car- 
rier concentration as evidenced by its larger 
conductance. The sample also shows little 
dependence on frequency of the real part of 
the conductivity and a small imaginary con- 
ductivity in this range. As the conductance 
is expected to drop off dramatically around 
the frequency of the scattering rate within 
the Drude model, this sample has a scat- 
tering rate greater than 1.0 THz. This is 



consistent with previous workpSEH Q n CVD 
grown graphene and on few-layer epitaxial 
graphene^!. 

IV. CONCLUSION 

The two different scattering rates found in 
our microwave and TDTS measurement show 
either the limitations of the Drude model for 
describing the fine details of electron trans- 
port in graphene at low frequencies, or alter- 
natively the inherent difficulties of extract- 
ing out the precise complex impedance of 
graphene. In the first case it may be that 
the details of scattering Dirac electrons give 
a conductivity lineshape that is not precisely 
Lorentzian. Then the slope of GijG\ can not 
be taken as a measure of r. In the second 
case the scattering rate may be underesti- 
mated by imprecisely removing the effects of 
the substrate contribution. It may be that 
this overestimates the size of the imaginary 
contribution to the impedance and thereby 
giving a larger contribution to G 2 - Future 
work will attempt to measure changes in 
the impedance by back gating the sample. 
It is likely that it will be possible to mea- 
sure changes in the impedance very precisely 
as the bias is swept from positive to nega- 
tive. Moreover, by going to the high bias 
regime/large conductivity regime we should 
be insensitive to influence of the substrate 



impedance and will be able to isolate the 
graphene conductance more precisely. 
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